International Journal of Trend in Scientific Research and Development (IJTSRD) 

Volume 4 Issue 6, September-October 2020 Available Online: www.ijtsrd.com e-ISSN: 2456 - 6470 

Phonon Modes and Elastic Properties of 
RE55A125Co20 Bulk Metallic Glasses 

P. H. Suthar 

Department of Physics, C U Shah Science College, Ashram Road, Ahmedabad, Gujarat, India 


ABSTRACT 


The present paper deals with the theoretical study of phonon modes and 
elastic properties of RE55A125Co20 (RE (Rare Earth) =Y, La, Ce, Pr, Nd, Sm, 
Gd, Tb, Dy, Ho, Er, Tm, Lu) bulk metallic glasses. The Hubbard-Beeby 
approach combine with Jani et al. local model potential has been utilized to 
compute the collective modes for present glasses. Three different local field 
correlation functions of Hartree (H), Taylor (T) and Farid et al (F) are 
employed to find the influence of the screening effects on the collective 
modes of RE55A125Co20 bulk metallic glasses. The calculated phonon 
dispersion curves show all broad features of the collective excitation. It 
observed that out of three local field correction function, Taylor and Farid 
et al screening functions give better agreement with experimental values 
compared to other screening functions. Overall, the present results of 
elastic properties are found to be good agreement with experimental 
results. 
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INTRODUCTION 

The bulk metallic glass (BMG) alloys have been considered 
as a special interest in the material science research over 
the last six decades. From the literature survey, bulk 
metallic glasses show unique structural features and 
outstanding mechanical, many novels, applying chemical 
and physical properties [1-9] because of glassy states, is 
associated widely with daily life, industry, materials 
preparation, organism preservation and a lot of natural 
phenomena [1-13]. 

Here, it has focused on rare earth base bulk metallic 
glasses (RE-Al-Co, {RE (rare earth) = Y, La, Ce, PR, Nd, Sm 
Gd, Tb, Dy, Ho, Er, Tm, Lu (lanthanide series)}) [3,4]. These 
glasses have distinguished itself from other metallic 
glasses due to their better glasses forming ability. 
Furthermore, rare-earth based bulk metallic glasses have 
concerned more attention for their unique physical 
properties. Such as mechanical properties, magnetic 
properties, elastic properties, thermodynamics, 
thermoplastic behaviour and bond properties of materials 
are one of the main topics for both scientific significance 
and practical application in condensed matter physics, 
materials science and engineering [1-13]. 

In regards of these Wang et al [3, 4] have reported the 
development of a series of the rare-earth (RE-Al-Cu) based 
bulk metallic glasses with desirable properties based on 
the correlation. The elastic properties of metallic glass can 
be roughly predicted to guide the selection of alloying 
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elements for developing a rare-earth base bulk metallic 
glass (RE-BMG) with higher glass forming ability (GFA). In 
this article. It is attention on the computation is carried 
out for phonon dispersion curve (PDC) and elastic 
properties likewise longitudinal sound velocity 'Vl', 
transverse sound velocity 'Vt , bulk modulus Bt, Poisson's 
ratio v, modulus of rigidity G, Young's modulus Y and 
Debye temperature 0d with the help of Hubbard Beeby 
[11] approach and model potential given by Jani et al [12- 
13] along with three different screening functions due to 
Hartree (H) [14], Taylor (T) [15] and Farid et al (F) [16] 
for the RE 55 AI 25 C 020 rare-earth bulk metallic glasses [3,4]. 

COMPUTATIONAL METHOD 

Hubbard and Beeby (HB) [11], Takeno and Coda (TG) [17, 
18] and Bhatia and Singh (BS) are main theoretical 
approaches used to calculate the phonon dispersion curve 
(PDC). One is the phenomenological theory of HB [11] in 
the random phase approximation, TG is the quasi¬ 
crystalline approximation technique with interatomic pair 
potential. And the third is by evaluation of force constants 
as was done by BS [19]. 

Out of this three theoretical approaches. The HB 
theoretical model has been utilized here to generate the 
phonon dispersion curves (PDC) of bulk metallic glasses of 
specific concentration (RE55AI25C020) in the present 
computation. 
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A. Model potential and Pair Potential 

The selection of model potential is an important aspect as 
it will describe the motion of valance electrons in a metal. 
For the present study, we have used single parametric 
local model pseudopotential of the following form in the q- 
space. [10,12,13]. 

^Oeffl 

exp(-l)(qrce//) 




cos{qrceff) - 
{sm{qrceff) + qrceffCOs{qrceff)} 


( 1 ) 


Here q, e and r^eff Is the valency, atomic 

volume, wave vector, charge of the electron and the 
parameter of the potential respectively. 


The effective interaction pair potential for the bulk 
metallic glass as a one metallic fluid, i.e, the concepts of 
effective atom [10], can be written as [10,13] 




^2 A 




O 


Oeff 




sin( qr) 


dq 


( 2 ) 


Here ^ and <-> are the effective valence and atomic 

eff ^'^Oeff 

volume of present system respectively. 

A quality which is equally important as the pair potential 
while studying a disorder system is pair correlation 
functions , which is computed theoretically from the 

effective pair potential. 

The energy wave number characteristics appearing in the 
equation (2) are written as [10-12] 

16 ^ J {l + [fi-# (q) - ill - f,ff C^)]} 




Here, iq) is effective bare ion potential it is given in 
equation (1), (<7) is the Hartree dielectric response 

function [12] and (^) is the local field correction 

functions to introduce the exchange and correction effects. 
The three exchange and correlation function due to 
Hartree (H) [without exchange and correlation effect) 
[15], Taylor (T) [12], and Farid et al (F) [13] are employed 
for the investigation of the relative effect of the local field 
correction function on PDC and elastic properties of 
RE55AI25C020 Bulk metallic glasses. 

B. Phonon Dispersion Curve and Elastic Constants 

The In the Hubbard and Beeby [13] approach, the 
longitudinal and transverse phonon frequencies are 
computed using [10,11,13]; 


ol{d)=^ 




3sm(^cr) 6cos(^cr) 6sin(^cr) 

qcr {qcrf {qo-Y 

3 cos(^cr) 6 sm( qcr^ 
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Here, co^ represents the maximum phonon frequency and 
is given as 
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and 


^(r) = exp 


I j 


C8J 


Here, p, Meff, g(r), (loeff, and Peff (q) be the number 
density, atomic mass, pair correlation function, atomic 
volume and energy wave number characteristic, 
respectively. 

In the long-wavelength limit, the phonon dispersion curve 
shows an elastic behaviour. Hence, the longitudinal vi and 
transverse vt sound velocities are also calculated by 
[ 10 , 11 ] 


V; = co^ 


vio; ^ ^1^10 


Various elastic properties are then determined by the 
longitudinal and transverse phonon frequencies. 

The bulk modulus Bt, Poisson's ratio modulus of rigidity 
G, Young's modulus Y and the Debye temperature 0d are 
calculated using the expression below [10,11] 


Bt = P\ 


G = pv] 


1- 2(v//v/) 

2- 2(vf/v/) 


Y = 2G(^ +1) 


and 
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Where, p is the isotropic density of the solid, 'h' is plank 
constant and ke is the Boltzmann constant. 

RESULTS AND DISCUSSION 

TABLE L: Input parameters and constants for 


Parameters -> 
BMG 1 

Zeff 

Doeff 

[au]3 

Rceff 

[au] 

Meff 

[amul 

Y55AI25CU2O 


161.29 

1.620 

54.09 

La55Al25CU20 


183.28 

1.607 

94.93 

Ce55Al25CU20 


184.66 

1.697 

95.59 

Pr55Al25CU20 


168.79 

1.712 

96.03 

Nd55Al25CU20 


166.50 

1.727 

97.86 

Sm55Al25CU20 


172.23 

1.730 

101.2 

Gd55Al25CU20 

2.085 

155.76 

1.675 

102.1 

Tb55Al25CU20 


158.47 

1.694 

101.9 

Dy55Al25CU20 


159.89 

1.868 

107.9 

H055A125CU20 


155.47 

1.694 

109.2 

Er55Al25CU20 


151.78 

1.730 

109.2 

Tm55Al25CU20 


150.88 

1.618 

110.5 

LU55A125CU20 


147.86 

1.651 

114.7 


With 
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All the calculations have been performed for the 
RE 55 AI 25 C 020 bulk metallic glasses. The input parameters 
and constants used in the present investigation of the 
phonon dispersion curve and elastic properties are 
tabulated in Table 1. Jani et al [14] model potential is used 
along with three different local field correction functions 
for the generate pair potential for the RE 55 AI 25 C 020 rare- 
earth bulk metallic glass system. From Fig. 1, it is observed 
that the position of the first minima is affected by the type 
of screening functions. The maximum depth in the pair 
potential in the present study is obtained due to F 
screening function for Ce 55 Al 25 Co 2 o BMG compared to 
other rare earth bulk metallic glasses. It is also noticed 
that the first zero for V{r = rQ)due to Hatree (H) Tq = 

13.76au~15.21au, for Taylor (T) [16] Tq = 

13.69au~15.08 au and for Farid et (F) [17]ro = 

13.71au~15.10 au are occurring for local field correction 
functions for RE 55 AI 25 C 020 rare-earth bulk metallic glasses. 
The oscillatory nature is also seen even at large r-region. 



Fig.l Screening dependence effective pair potential of 
RE55AI25C020 BMG 


The computed results of the longitudinal phonon 
dispersion curve [o)l) and transverse phonon dispersion 
curve [ o ) t ) generated using HB scheme along with three 
different screening functions are shown in Fig. 2 and Fig. 3 
for the RE55AI25C020 rare-earth bulk metallic glasses 
respectively. In an absence of experimental data for 
structure factors, we can offer reverse comment that the 
structure factor S (q) will show it's at these qs. It is 
understandable from the figures that the nature of peak 
positions are not much affected by different screening 
functions, but both the transverse and longitudinal 
frequencies show deviation by using T and F with respect 
to H screening functions. It is observed that the present 
fallout of phonon modes with correlation function due to T 
function is positioned between those due to H and F 
screening functions. The highest value of first peak point 
and lowest values of first peak point are seen for 
Y55AI25C020 and Ce55Al25Co2o BMG respectively, and rest of 


rare earth bulk metallic glasses peak point fallout between 
these two bulk metallic glasses. It is clear from the figures 
that the oscillations are prominent in the longitudinal 
mode as compared to transverse mode, which indicates 
that the collective excitations at larger wave vector 
transfers due to the dispersion of longitudinal excitation. 
On the other hand the transverse modes undergo larger 
thermal modulation due to the anharmonicity of the 
vibrations in the bulk metallic glasses. In the long 
wavelength limit the dispersion curves are linear and 
confirming characteristics of elastic waves. 



q^) 

Fig.2 Screening dependence longitudinal collective 
modes of RE55AI25C020 BMG 



q(A-^) 

Fig.3 Screening dependence transverses collective 
modes of RE55AI25C020 BMG 


Presently computed elastic properties for RE 55 AI 25 C 020 
rare-earth bulk metallic glasses are presented in Fig. 4 to 
Fig 10 for longitudinal sound velocity, transverse sound 
velocity, bulk modulus Bt, Poisson's ratio v, modulus of 
rigidity G, Young's modulus Y and Debye temperature 0d 
along with experimental values respectively. The 
experimental values computed using the ultrasonic 
method is reported by W. H. Wang [3]. From these figures, 
one can see that the values are obtained due to T [16] and 
F [17] screening functions are very close to one another as 
compared to the H screening function. From the Fig 4 and 
Fig 5, the comparison shows that presently computed 
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results are obtain lower values for longitudinal velocity 
but better agreement with transverse sound velocity 
respectively. From Fig 6 , one can see that outcome of bulk 
modulus show the underestimate computed using present 
model potential. But the F screening function shows the 
less deviation compared to H and T screening functions. 
From Fig. 5, Fig. 7, Fig.8 and Fig. 10 are shown the excellent 
agreement with available experimental values for 
transverse sound velocity, modulus of rigidity 'G', Young 
modulus T and Debye temperature 0d respectively. The 
experimental data are not available for Y 55 AI 25 C 020 , 
Ce55Al25CO20, Nd55Al25CO20, Sm55Al25CO20 and Gd55Al25CO20 
bulk metallic glasses. So, at present we do not put any 
concrete remark for this bulk metallic glass. The computed 
elastic properties using T and F screening functions give 
excellent agreement particular for Tb 55 Al 25 Co 2 o, 
Dy55Al25Co2o and H 055 AI 25 C 020 bulk metallic glasses 
compare to in the series of other rare-earth bulk metallic 
glasses. From the Fig 4 to Fig 10 it is observed that the 
result generated due to HB-approach along with Jani et al 
model potential are found to be in acceptable agreement 
with the experimental results given in respective 
literatures. 



RE^AI^Co^BMG 


Fig.4 Longitudinal sound velocity (Vl) RE55AI25C020 
BMG along with available experimental values [3] 



RE,AI,Co,BMG 


Fig.5 Transverses sound velocity (Vt) RE55AI25C020 
BMG along with available experimental values [3] 



RWo^bmg 

Fig.6 Bulk modulus (Bt) of RE55AI25C020 BMG along 
with available experimental values [3] 



RE^Al Co.,BMG 

fib 25 2D 

Fig.7 Modulus of rigidity (G) of RE55AI25C020 BMG along 
with available experimental values [3] 



Fig. 8 Young's modulus (Y) of RE55AI25C020 BMG along 
with available experimental values [3] 
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RE,AI,,Co,,BMG 

Fig.9 Poison ratio (u) of RE55AI25C020 BMG along with 
available experimental values [3] 



RE55AI,,Co^BMG 


Fig.lO Debye temperature (0d) of RE55AI25C020 BMG 
along with available experimental values [3] 

TABLE 11 : RE55AI25CU20 percentile deviation with 


respect to experimental values 


Elastic 

Properties 

H[15] 

T[16] 

F[17] 

Longitudinal 
sound velocity 
'Vl' 

17.96%- 

22.34% 

2.33%- 

9.40% 

1.77%- 

8.50% 

Transverse 
sound velocity 
Vt' 

6.89%- 

15.79% 

0.14%- 

11.16% 

0.45%- 

11.77% 

Bulk modulus 
'Bt' 

38.35%- 

55.68% 

18.68%- 

30.12% 

18.56%- 

28.72% 

Modulus of 
rigidity 'G' 

17.18%- 

33.31% 

0.52%- 

23.62% 

1.14%- 

24.98% 

Young modulus 
T 

20.02%- 

36.99% 

0.69%- 

16.22% 

0.51%- 

17.74% 

Poison Ratio 'u' 

18.30%- 

23.54% 

18.3%- 

23.54% 

18.30%- 

23.54% 

Debye 

temperature 0d 

9.10%- 

33.48% 

0.29%- 

16.50% 

0.24%- 

15.63% 


The maximum and minimum percentile deviation of the 
various local field correlation functions with respect to the 
available experimental data is tabulated in Table 2. From 
this Table 2, the local field correction functions due to H 
(without exchange and correlation) give maximum 
percentile deviation compare to T and F functions and F 
functions gives the minimum deviation compare to other 
screening functions. So that we put the remark F screening 
function is very close to existing experimental data 
compared to H and T screening functions. 

However, the improvement in the present finding may be 
achieved either by incorporating other forms of exchange 
and correlation effect or by suggesting the modification in 
determining the parameters of the potential. From this 
calculation, it is observed that the presently computed 
data serve as a guideline for further research in this 
direction. 
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